Coherent output of photons from coupled superconducting transmission line 

resonators controlled by charge qubits 
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We study the coherent control of microwave photons propagating in a superconducting waveguide 
consisting of coupled transmission line resonators, each of which is connected to a tunable charge 
qubit. While these coupled line resonators form an artificial photonic crystal with an engineered 
photonic band structure, the charge qubits collectively behave as spin waves in the low excitation 
limit, which modify the band-gap structure to slow and stop the microwave propagation. The 
conceptual exploration here suggests an electromagnetically controlled quantum device based on 
the on-chip circuit QED for the coherent manipulation of photons, such as the dynamic creation of 
laser-like output from the waveguide by pumping the artificial atoms for population inversion. 

PACS numbers: 85.35.Ds, 73.23.Hk, 42.70.Qs, 03.67.-a 



I. INTRODUCTION 

Recent experiments with on-chip all optical setups [H, 
[3, S 0] have displayed slow light phenomenon similar to 
that due to the electromagnetically induced transparency 
(EIT) I, [|. Here, a physical mechanism was is pre- 
sented using a model with a coupled resonator optical 
waveguide, which behaves as a photonic crystal with a 
band-gap spectrum. The coupling of each resonator to 
some external cavities can shift the resonant spectral line 
and compress the bandwidth, and thus stop or store the 
propagating light pulses [J, 0, Q ■ 

Motivated by this progress, both in experimental and 
theoretical aspects, we propose and study a hybrid struc- 
ture with a cavity waveguide interacting with two-level 
artificial atoms. Here, the novel control mechanism 
for coherent transmission of microwave photons in the 
waveguide is to utilize the collective excitations of the 
atoms, which can be described as quasi-spin waves 
in the low excitation limit. As illustrated in Fig. [1] we 
suggest a co-planar on-chip setu p ba sed on the supercon- 
ducting circuit QED [1, i, 0, [111, [H [H, [H, [S [H]: 
each cavity in our proposed setup has been experimen- 
tally implemented as a superconducting transmission line 
resonator; the spatially distributed artificial atoms are 
the biased Cooper pair boxes (charge qubits), which 
play the same role as the external cavity for controlling 
light in the waveguide of the on-chip all optical experi- 
ment [Ei,Sll. 

Due to its engineered photonic band structure, such 
a cavity array, coupled to quasi-spin wave excitations, 
can result in much richer quantum coherent phenomena. 
We show that the quasi-spin waves of charge qubits can 
controUably affect the engineered photonic band struc- 
ture so that the modified dispersion relation results in 
a slow (and even zero) group velocity of photons prop- 
agating along the waveguide of coupled line resonators. 
The spin-wave excitation can also drive a coherent state 



of photons in the coupled resonator waveguide. Some 
pumping methods can make a population inversion to 
produce a laser-like output above the threshold. 

This paper is organized as follows. In Sec. [Tll we 
present our hybrid system: an array of coupled-line- 
resonators based on a superconducting circuit. Each res- 
onator is coupled to a biased Cooper pair box (CPB). The 
array of coupled-line-resonator exhibits a band structure. 
In Sec. mil in the large- and low-excitation limits, the 
hybrid system is modeled as two coupled boson models, 
with one for the photonic band and the other for the 
spin wave of N CPBs. Here, the slow light phenomenon 
is found by controlling the detuning and the coupling 
strength. The dynamic creation of the laser-like output 
from the array of coupled-line-resonator is proposed in 
Sec. IIVI and |Vl This is because two kind of bosons cou- 
ple linearly with each other and the population inversion 
of the artificial atoms is pumped. Sec. IVII presents our 
conclusions. 



II. THE COUPLED CAVITY QED BASED ON 
SUPERCONDUCTING CIRCUIT 



Now we consider a superconducting quantum cir- 
cuit including A^ CPBs and an array of couplcd-line- 
resonator. Recently, most experiments have demon- 
strated the possibility to fabricate a superconducting 
qubit array ^181, [l^ . As shown in Fig. [2 the array 
of coupled-line-resonators is constructed by cutting the 
superconductor transmission line into N segments. And 
the A^ CPBs are made of A^ dc SQUIDs (superconduct- 
ing quantum interference device) which consists of two 
tunnel junctions. 

With a proper biased voltage Vg, each CPB behaves 
as a two- level system (charge qubit) [l^l- Typically, the 
model Hamiltonian for a charge qubit can be written as 
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FIG. 1: (Color online) Configuration of the setup for control- 
ling light propagation in a coupled-line-resonator waveguide 
(a) by coupling to charge qubits (b). The coupled (line res- 
onators) - (charge qubits) system can behave similarly to the 
cavity QED for a single atom interacting with a single-mode 
cavity (c). 

in Ref. 

H = ^B,a, - ^B.,a^ (1) 

with the level spacing of the charge qubit 

= 4.Ec i2ng - 1) (2) 
and effective Josephson energy 

B.,^2Ejcos(^7T^y (3) 

The quasi-spin operators cr^ and ax are defined in the 
charge qubit basis (|0) and |1)), where and 1 repre- 
sent excess Cooper pairs on the superconducting island, 
respectively. Here, Ec and Ej represent the charge en- 
ergy and the Josephson energy, respectively. Through 
Ug = CgVg/2e, wc cau tune the degenerate point of the 
charge qubit by controlling the gate voltage Vg applied 
on the gate capacitance Cg. In addition, represents 
the magnetic fiux through the SUID loop induced by the 
externally applied magnetic field and <I>o = h/2e is the 
flux quantum. 

Based on the above results, we can write the Hamilto- 
nian including N CPBs as 

^A = ^y^(|e,)(e,|-|5,)(<?,l) (4) 
j 

where the eigen-frequency lua — yB^+B^. Based on 
recent experiments, the charge energy can be taken as 
Ec — 29.5 GHz, while the Josephson energy Ej = 8 
GHz. For simplicity, all qubits are assumed to be iden- 
tical and biased off the degenerate point. The energy 
eigenstates 

|e,>=cos(^0 |O,)-sin(^0 |1,) (5) 



and 

l5.)=sin(^0|O,)+cosQ)|l,), (6) 

are the superpositions of charge eigenstates \0j) and 
and 9 — arctan{Bx/ B^)- 

In our setup, each qubit is placed at the position of the 
anti-node of the standing wave field in each transmission 
line resonator. The London equation provides the van- 
ishing boundary conditions for the quantized electromag- 
netic field at the two ends of each line resonator. Thus, 
the quantized magnetic field vanishes at those anti-nodes 
and the qubits are only coupled to the electric compo- 
nent. Here, the gate voltage applied to the capacitance 
of jth Cooper pair box is given by — Vg + , where 
stands for the quantized part of the voltage. Let only 
the n = 2 mode be activated. Then the gate voltage is 
quantized as [H, [l^ 

= (7) 

at its maximum while the magnetic field vanishes. Here, 
we only assumed a single mode of the quantized electric 
field of frequency co with creation (annihilation) operator 
{dj). L is the length of the line resonators and C is 
the capacitance per unit length. 

There exists the coupling J between two neighbor line 
resonators through the dielectric material. Then the 
Hamiltonian of the coupled line resonators waveguide is 

He = ujY^ d]dj + JY^ {d]dj+i + a]+iaj) , (8) 
j j 

where J depends on the coupling mechanism, e.g., the 
magnetic penetration depth. 

Following recent experiments [l3| , we can take the fre- 
quency of the line resonators uj = 10.0 GHz, the length 
of line resonators L = l.O cm and the capacitance per 
unit length C — 0.13 fF //^m. We now make the rotat- 
ing wave approximation to write down the interaction 
Hamiltonian [H, [li| 

N 

Hi = gY,{d,\e,){g,\+R.c.), (9) 
between the qubits and the fields, where 




is the coupling strength between the resonator and the 
artificial two-level atom, and Cs is the sum of the gate 
capacitance and the capacitance of the tunnel junction. 

In practical experiments, the coupling constants g and 
J should depend on the position of the qubit. For sim- 
plicity, in this paper we take a uniform g and J. The- 
oretically, the fluctuations of the coupling constants are 
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assumed to be innocuous and do not change the results of 
this paper, quahtatively. The above model has been used 
to demonstrate the photon blocked effect [H, [l^, \2^ . 
which leads to the Mott insulating effect for the polari- 
tons formed by dressing atoms with a cavity field. 

Although our proposed configuration setup is an ex- 
tension of the single cavity proposed in [l^, [H, 
to N coupled cavities, there are significant difference be- 
tween them. In experiments [l3| , to demonstrate the 
typical cavity QED character of a superconducting quan- 
tum circuit for quantum computing, only a single cavity 
is used rather than a coupled cavity array. A single cav- 
ity possesses discrete photon modes with large frequency 
spacings, but an array of coupled- line-resonator possesses 
a band-gap spectrum and can transmit a wave packet of 
light. A single cavity coupled to many charge qubits in a 
similar superconducting circuit has been modeled [1^ to 
probe the dynamic behavior of quantum phase transition 
j20 |. Obviously, different from the setup in [13, 14], our 
setup uses many coupled transmission line resonators to 
realize an electromagnetically controlled quantum device 
for coherent control of photon transmission. 



III. TWO-MODE BOSON MODEL FOR THE 
SPIN- WAVE DRESSED PHOTONIC BAND 



Here, we have introduced the polariton operators 0] 

Pfc = cosOkak + smOkhk; (14) 
Qk = sin^fcOfc - cos6kbk. (15) 

which are linear combinations of the quantized electro- 
magnetic field operators and atomic collective excitation 
operators of quasi-spin waves. The mixing angle 9 is de- 
termined by 

tan9k = 2g/{nk~ujA)- (16) 
The total excitation number for the kth mode is given by 



Nk = PlPk + QlQk = alak + blk 



At; 



(17) 



The dispersion relation for the photonic band is obtained 
as 



where 



e±k = ^Dk ± Efc, 



^Dk = ^(^fc +^a), 
Sk = ^y/ink-LUA)^ 



(18) 



(19) 



Now we consider the low-excitation limit that a few 
charge qubits are populated in their excited state. The 
crucial issue here is to use the collective operators |7i] 



N 



(11) 



and its conjugate Bk — (^1)'' to describe the spin wave 
excitation of the charge qubit array, where k = 2'Kn/£N 
with n = 0, 1, — 1. In the large A^ limit with low 
excitations, i.e., {J2k -^l-^k) A^, these collective exci- 
tations behave as bosons, since the usual bosonic commu- 



Bk.Bi 



Skk' can be approached when 



tation relation 
A^ ^ oo. 

The Fourier transformation ak — X]j Gxp{ik£j)aj / ^/N 
diagonalizes the coupled resonator Hamiltonian as 



k 

to give a dispersion relation with band structure 



fifc = w + 2Jcoa{k£), 



(12) 



where £ is the site distance. Then the normal modes 
of the hybrid system with Hamiltonian H — J^k Bk are 
characterized by ak and bk = lim^r^oo Bk- The evolution 
of each mode is governed by the following Hamiltonian 

Hk = ^ka\ak + g (akbl + H.c.^ + ujAblbk (13) 
= riDkNk + ekiPlPk - QlQk). 



Since A^fc commutes with Hk, the number of excitations 
A'^fc is conserved, while the number of different type exci- 
tations a^flfc and bjjjk are mutually convertible by adjust- 
ing the coupling strength g and the detuning S — uj — oja 
between the artificial atom and the resonator. Due to 
the coupling between artificial atoms and the resonators, 
the origin band structure for the array of coupled-line- 
resonators is splitted into two in the one excitation sub- 
space. Obviously there exists a gap between two bands 
for nonvanishing g. 

Notice that the band structure (fT9)) is only available in 



the low excitation Hmit {J2k ^k'^k) ^ A^ ■ In the following, 
we focus on the single excitation subspace. For S > 0, the 
corresponding dressed spectrum e±k — fiufe ± shows 
a bandwidth narrowing effect on the coupled line res- 
onators due to its couplings to the charge qubits. This 
is because the bandwidth W- — \e-k=o — e-fc=7r| of the 
low-band changes from 2 J to 2J — A, where 



A = A{J)-A{-J) 



and 



A{J) 



jy + g-'. 



(20) 



(21) 



Meanwhile, the bandwidth W+ = |e+fe=o — e+fc=7r| of the 
upper band changes from 2 J to 2J + A. Obviously when 
S = 0, the bandwidths and W+ of the two bands are 
equal to 2 J; when 6 > {S < 0), W-^ (W^+) is narrower 
than 2 J while W_|_ {W-) is wider than 2 J. Then the wave 
packet of light in the lower band can be adjusted by S. 
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The couplings also shift the central spectral line from lu 
and loa to e±fc=^/2 respectively. This just recovers the 
same result obtained in Ref. [J [3, S 13 • 

Then we can obtain the group velocities v±{k) = 
de±k/dk for the lower and the upper bands as 



= Re < Jls\n{M) 



1± 



A + 2Jcos(H) 



2Jcos(H)]2 +4g2 

"(22) 

where A = 5 + in. Here, n ~ rj — ^ \s the difference 
between the damping constant 77 of each cavity and the 
decay rate 7 of the qubit, where rj is caused by the cav- 
ity loss and 7 = Tj~^ + T^^: in experiments, the en- 
ergy relaxation time Ti of a superconducting qubit is 0.1 
to a few microseconds, and the dephasing time T2 be- 
tween the ground and the excited state is a few dozens 
of nanoseconds. It can be seen that the group velocity 
becomes independent of loss when 77 = 7. At the band 
center k = n/{2t), the group velocity for the lower band 
becomes 



= Re <^ 1 



A 



v/A2 + 4.g2 



(23) 



In the large detuning case, i.e., 5 ^ 2|g|, the group ve- 
locity reaches its minimum 



t;_(-7r/(2^)) =Re [Jig^/K^ 







(24) 



corresponding to the lower bandwidth compressed. In 
the case of 6 ^ — 2|g|, the group velocity for the lower 
photonic band reaches its maximum 



V-{-it/{21)) w 2Jt. 



(25) 



And the lower band has a large bandwidth, which ac- 
commodates the entire pulse bandwidth. In Fig. [21^ a) we 
have plotted the group velocity of the lowest band as a 
function of the level spacing uja of the two-level artificial 
atom and the coupling strength g dX k = — 7r/(2€). It can 
be seen that when a large detuning occurs, the group ve- 
locity is minimum at large blue detuning and is maximum 
at large red detuning. Hence, for a microwave pulse as 
a superposition of many fc-states, its distribution in the 
fc-space can be entirely contained in the energy band by 
setting 5 <^ — Slg]. Therefore the microwave pulse can 
be stopped by adiabatically tuning the detuning from 
5<^-2\g\ t0(S>2|5|. 

In the near- resonance case 5 ~ with strong couplings, 
the group velocity in the lower band 



w_ = Re<^ J£sin(H) 



Jcos{k£) + IK 



is reduced to zero approximately within the range of 
k satisfying J cos{k£) ^ g. For mode k, which satis- 
fies J cos{k£) <C —g, the corresponding group velocity 
becomes V- « 2J^sin(fc^). For mode k, which satis- 
fies — g <C Jcos{k£) ^ g, the group velocity becomes 




FIG. 2: (Color online) (a) the group velocity Vg of the 
lowest band as a function of uja and coupling strength g 
at k — ■k/(2£) for — y, and u = J. It can be seen 
that when 5 ^ 2g, the group velocity reaches its minimum 



Vg{~n/{2£)) fa 0; when S < 
its maximum Vg{—n/{2£)) « 



-2g, the group velocity reaches 
2Jl. (b) the group velocity of 



the lowest band as a function of k and coupling strength g in 
the case of resonance. Here, g and uja are in units of J, and 
k is in units of the lattice constant £. 



V- ~ J^sin(fc^), which are illustrated in Fig. [2l^b). It 
tells us that when a microwave pulse inputs into such a 
coupled-line-resonator waveguide, some components will 
be stopped completely by adjusting the coupling strength 
g, while others still pass through. So the microwave pulse 
is distorted. Hence, in the case of resonance, one cannot 
obtain the whole information that the wave packet car- 
ries. 



IV. COHERENT OUTPUT OF SLOW LIGHT 

The above boson model describes the linear couplings 
between two kinds of boson modes. If one can prepare the 
state of the charge qubit array with population inversion 
through coherent pumping, the dynamic evolution will 
drive the coupled-line-resonator mode to output a laser- 
like light in a coherent state. Such pumping mechanism 
with superconducting qubits has been explored most re- 
cently for superconducting flux qubits [23]. 

Now we assume the charge qubits are prepared co- 
herently in excited states. Then the initial state of the 
spatially-distributed atomic ensemble is a fc-mode coher- 
ent state 

\ak) = D{ak)\G) = D{ak)\gi) ® I52) «) . . . ® Isw) (26) 

where the displacement operator 

D{ak)^exp{akbl-albk), (27) 

describes a superposition of rt-quasiparticle excitation 
states \nk) of mode-fc. Here, the one quasiparticle excita- 
tion characterizes the spatially-distributed qubit popula- 
tions with definite phases, a quasi-spin wave in the charge 
qubit array. Let the total system initially start with a co- 
herent state of atomic ensemble and the vacuum state of 
the coupled-hne-resonator array, i.e. |'(/'(0)) = |0) <S) \ak)- 
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Here, |0) means no photon contained in the coupled-hne- 
resonator array. After time t, the initial state is evolved 

into m)) = U{t)\i;m. 

In order to find an explicit expression for \ip{t)), we 
formally rewrite the quantum state at time t as 

m)) = U (t) D{ak)U-^ {t) U {t) \G) |0) (28) 

where U (t) = exp(ii iJ^) is the time evolution op- 
erator. Because the number of excitations in the total 
system is conserved and 10)10) is the ground state of the 
system corresponding to the zero eigenvalue, the quan- 
tum state defined by Eq. (^5]) becomes 



(29) 



\m) ^Uit)D{ak)U-'{t)\G)\0), 



which is completely determined by the time dependent 
displacement operator 



D{t) = U (t) D{ak)U-^ (t) = exp A {-t) 



(30) 



where A {-t) = U (t) AU'^ {t). 

With respect to the polariton operators P and Q, 
the initial coherent state \ak) of the fcth mode can be 
rewritten based on the displacement operator D{ak) = 



exp 



A{0) 



where 



A{0) = akPl sin Ok - akQl cos 0k - h.c. (31) 

is an anti-Hermitian operator. Since polaritons are the 
eigenexcitation of the total system, their creation opera- 
tors PI and qI. evolve according to the eigenfrequencies, 
i.e., 

U{t)PlU-^{t) = ple-^(^o,+e,)t ^ 
U{t)Qlu-^{t) = Qle-'^^""'"-^'. 

Accordingly, the operator A{t) can be explicitly obtained 
as 



A{t) 



ptg-«(0^.+e.)igj 



+alQke 



i(0£ 



sm t/fc 
' 008 6*^ — h.c. 



(32) 



Transformed back to the original representation with op- 
erators Obk and 6fc, the displacement operator D {t) be- 
comes the product of two displacement operators, i.e.. 



D{t)^D [ak m D [Pk (t)] , 
where the factor 

D [ak {t)] = exp \ak (t) bl - al (t) bk 
acts on the atomic excitation state while 

D[Pk{t)]^e^p \f3kit)ai-f3l{t)ak 



(33) 
(34) 

(35) 



acts on the state space of the coupled-line-resonator ar- 
ray. 



Therefore, the state \il^{t)) can be factorized as \ip{t)) 
\Pk{t)) ^ |Q!fc(^)) with the time-dependent amplitudes 



-''"'sin^ 



ak{t) = uke-''^""' [e'^"' cos" Ok + e 
(3k {t) = -te-'""'''aksm{ekt)sm{29k), 



0fcO3,6a) 
(36b) 



Here the coherent state \Pk{t)) localizes in the fc-mode, 
and it actually is a spatially multi-mode coherent state; 
the periodically modulated amplitudes Pj{t) mean a 
quasi-classical wave packet of photons with distribution 
Pj = \/3j{t)\'^. The above argument means that the ini- 
tial coherent input of the atomic excitation can result in 
a coherent output in the photonic mode, which is char- 
acterized by a coherent state \Pk{t)). This just displays a 
laser-like behavior for the photons output from the cou- 
pled line-resonator waveguide. 



V. LASER LIKE PROCESS WITH 
POPULATION INVERSION 

The above intuitive discussion shows an obvious laser 
behavior, but we need the population inversion imple- 
mented by some coherent pumping. We also need to con- 
sider the threshold condition for the lasing in the coupled 
line-resonator waveguide. To this end we study the co- 
herent radiation in the coupled line-resonator waveguide 
stimulated by the the artificial atoms with a collective 
coherent excitation. 

Ignoring the fluctuations due to the couplings to the 
thermal bath, the dynamic variable of this system obeys 
the following equations 



dt&k 
dtBk 

dtSo 
where 



~i{flk - ir])ak - igBk, (37a) 
-i{tJA~h)Bk+ijj:'^Sk'-kak', (37b) 

k' 

r (^Nda - So) ~ i2gY^ (plak ~ a+Bd^7c) 

k 



N-1 



Sk'-k = 2^ e I )J<Tj 



j=0 



and So is defined hy k' = k. By neglecting the scattering 
from k to k' due to lar ge N , one can write down a system 
of laser-like equations [2^ 



dtak = -ii^lk - «'7)afc ~ «.9^fe: 
dtbk = -i{ujA ~ n)bk + igncik/N, 
dtfi = r {Ndo -~n)~ i2gY, - a+bk ) , 

k 



(38) 



where do is the input rate for equilibrium inversion. 
The term T {Ndo — n) in Eq. ([55)) is phenomenologi- 
cally introduced to characterize the role of pumping 
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by some population inversion for the excitation number 

3 to equilibrium, 
time-dependence in 



fi = ^bj,bf^; T is the relaxation time to equilibrium 



Next, we eliminate the fast 
Eq. pS]) by the substitutions 



and 



bk = bk exp(-irifct) 



(39) 



(40) 



We lock the cavity mode in the coupled-line-resonator in 
a resonance frequency Q.k ~ ^A, that is, exp[i(w^ — Q,k)] 
varies slowly. When the relaxation time of charge-qubits 
is much smaller than the relaxation time of the cavity as 
well as that of the population inversion, i.e., 7 3> T ^ 77, 
we can adiabatically eliminate bk and n by setting the 
corresponding time derivatives to zero in the obtained 

equations about ak, bk and fi from Eq. ()38|) . Finally we 
obtain the equation of motion for hk 



dtCLk = [dog^L - ri)ak 



ivr 



where the Lorentz distribution 



1 



7 + i{uJA - ^k) 



(41) 



(42) 



shapes the spectra of coherent output. It is exactly the 
laser equation [28] defining a threshold do — rj/{g^L) for 
the laser-like output in the coupled-line-resonator waveg- 
uide by the coherent pumping with an input rate do- 

This dynamical lasing behavior can be described by 
the non-linear equation 



(43) 



about the order parameters x — {a\ak\a). Here, \a) is a 
coherent state with real number a\ the coefficient of the 
linear term of Eq. (|43|) 



(44) 



a = rj — dog L 



describes the amplification effect of the light field when 
cr < 0; the coefficient 



A = 



A^dpg^L \L\ 
NT 



(45) 



represent the nonlinearity of the effective theory obtained 
by averaging the atomic excitations. The nonlinear coef- 
ficient competes with the parameter a to realize a lasing 
"phase transition" . As illustrated in Fig. [31 the solutions 
of the non-linear equation obviously possess a critical be- 
havior near the threshold ct = . The solution of Eq. ((43)) 
gives 



[cexp(2cri) - A/cr]"^ 
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FIG. 3: (Color online) The probability amplitude |a;p of out- 
put photon as a function of a and time. When a < 0, there is 
a laser-like output depicted by the stable solution for t ^ 00. 
t is in units of seconds "s" and a is in units of s~^. 



for a ^ 0; and 



\x\^ = (2At + c)' 



(47) 



for (7 = 0, where c is a constant determined by the initial 
state. The stable solutions that 



-cr/A 



(48) 



(46) 



for cr < (above the threshold) and |a;| = for cr > 
(below the threshold) mean a laser-like output of the 
coupled-line-resonators. Since the expectation value of 
(ik does not vanish, we can understand the above laser- 
like effect as a pumping-induced symmetry-breaking. 

Together with the intuitive argument in the last sec- 
tion, the analysis in this section definitely shows the exis- 
tence of a threshold for laser-like behavior in the present 
artificial system. However, to reduce the threshold by 
overcoming decoherence, including dissipation and de- 
phasing, is still a challenge to implement a laser of slow 
light in practical hybrid systems. 



VI. CONCLUDING REMARKS 

In conclusion, we have conceptually proposed an elec- 
tromagnetically controllable quantum device based on 
superconducting circuit QED for the coherent manipula- 
tion of photons. It can realize a laser-like slow microwave 
output from a coupled-line-resonator waveguide by con- 
trolling each cavity connected to a charge qubit. Our 
studies are motivated by the all-optical experiment for 
stopping light in a coupled resonator waveguide, which 
is coupled to other cavities, but here we show that some 
results for the coupled resonator waveguide can be ob- 
tained by replacing the coupling cavities with a more 
practical system, a spatially-distributed array of charge 
qubits. Finally, we would like to mention that, without 
the limitation of the superconducting system, the present 
analysis represents a universal setup for the coherent ma- 
nipulation for light or microwave propagations in some 
all-optical or electromagnetic-optical system. 
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